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Abstract 

We study a class of interacting, harmonically trapped boson systems at angu- 
lar momentum L. The Hamiltonian leaves a L-dimensional subspace invari- 
ant, and this permits an explicit solution of several eigenstates and energies 
for a wide class of two-body interactions. 
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The study of vortices in Bose-Einstein condensates of trapped atomic vapors is of con- 
siderable experimental [||,|| and theoretical @-13] interest. Present experiments and much 



theoretical work focus on the Thomas- Fermi regime of short coherence length The 
opposite limit of perturbatively weak contact interactions is not yet realized experimentally 
but is of theoretical interest as well p|-p!3|. The case of attractive interactions was studied by 
Wilkin et al. . Mottelson developed a theory for repulsive interactions fl^ , and Kavoulakis 



et al. [O] and Jackson et al. [O] compared mean field and exact numerical results. Exact 



diagonalization techniques by us showed that the ground state energy depends linearly on 
angular momentum and led to an analytical expression for the ground state wave function 
Tl| . Recently, the form of this wave function and its eigenenergy were confirmed in analytic 



calculations |[15Hl7|. However, we still lack insight about why certain eigenstates turn out 
to be simple analytical functions. 

We show in this paper that there is a subspace structure that explains these findings. 
Interestingly, this structure is not limited to the two-body contact interaction considered in 
previous investigations but present in a rather large class of Hamiltonians. Let us consider a 
system of harmonically trapped bosons in two spatial dimensions Q at angular momentum 
L. The trap Hamiltonian is 

H-\T.{r^] + r]-2), (1) 

i=i 

where r| = x^+i/j. Let us use complex single-particle coordinates zj = Xj+iyj, j = 1 . . . N 
and fix the angular momentum to be L. The eigenfunctions of the trap Hamiltonian are of 
the form 



i){zi, . . . , zn) = (pizi, . . . , zn) exp --HkiH (2) 




and have degenerate energy E = L. In eq.(|^) 4>{zi, . . . , zn) denotes a homogeneous polyno- 
mial of degree L that is totally symmetric under permutation of particle indices. Suitable 
basis functions for such polynomials are products 

4>{Zi, ...,Zn)= CAiCAa • • • 

where ex denotes the elementary symmetric polynomial 

^xi^ly ■ ■ ■ 1 ^n) — ^pi^p2 ■ ■ ■ ^p\-i 

l<pi<P2<...<px<N 

and {Al, A2, . . . , Afc} is a partition of L into at most integers, e.g. Aj = L with 

Al > A2 > . . . > Afc > and k < N. We set cq = 1. From here on we omit the notation 
of the ubiquitous exponentials in the wave function and restrict ourselves to the regime 
< L < N. In this regime, N is simply a parameter of the problem, and the dimension of 
Hilbert space equals the number of partitions of L into integers. 



^The three dimensional problem is essentially two-dimensional, see e.g. ref. IT 
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The operator for the total angular momentum is 



N 

L = T.^A^ (3) 

where dj denotes the derivative with respect to Zj and acts only onto the polynomial part 
of the wave function. Another important observable is the angular momentum of the center 
of mass. Its operator is given by 

= z^ D^, (4) 
where Zc = Ci/N denotes the center of mass and 

TV 

Dc = T.dj. (5) 

Again, it is understood that derivatives act on the polynomial part of the wave function only. 
Lc has eigenvalues 0,1,2,...,L — 2,L and commutes with the harmonic trap Hamiltonian 
(H) and the angular momentum (^). We further have the following commutator relations 

[L, Dc] = -De, [L, Zc] = Zc, [Lc, Dc] = -De, [Lc, Zc] = Zc. (6) 

The introduction of perturbatively weak interactions lifts the degeneracy of the trap 
Hamiltonian. In what follows we are interested in two-body interactions of the form 

L 

^ = 9 ^ Cm^mj 
m=0 

A^= E {z, - z,r {d, - d,r , (7) 

l<j<jr<Af 

where the derivatives act only onto the polynomial part of the wave function. We choose 
gN <^ 1 to be in the perturbative regime. Note that a large class of two-body interactions 
is of the form ([^. Consider for instance the two-body potential V{r) = r"^^. Comparison of 
matrix elements shows that 



H^r^-r^r) =lj:[{zl-z*^){z,-z,) 



3) ^3) 



\m=0 

Thus, any potential that depends analytically on the squared inter-particle distance can 
be written in the form of eq. (0) . Further examples include zero range potentials of the form 

V{r) = 27vg \ao6{r) + aiV^(5(r) + a2V^6{r] 



3 



For Cm = (— l/2)"*/m! the interaction (|^) corresponds to the two-body contact interaction 
V = 2'Kg5{r). This can be seen by computing and comparing matrix elements of these two 
different representations of the interaction. It is also instructive to analyze the action of 
the operator (^) on the wave function. Inserting the appropriate coefficients for the contact 
interaction into the Hamiltonian (0) we obtain a Taylor series. Thus, 

^ 1 / Z ■ — Zi\"^ m 

—] ( ^ 9 ) ~ ^jT 0(^1, . . . , ^i, . . . , ^j, . . . , ^iv) 

m=0 \ ^ / 

L ( '\\n ~ \ m+n 

= E bir drd-^{z,,...,z,,...,z,,...,z^) 

= (p \Zi, . . . , Zi , . . . ,Zj , . . . ,Zn 
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^1 5 • • • 5 ~ ) • • • ) ^^t; , ■ ■ ■ ,zn 



~\~ Zj ~\~ Zi 



In ref. |]T^ we presented an alternative differential operator with this effect on wave functions. 

The operators Am commute with L, Lc, Zc and D^. Eq.(@) thus implies that the applica- 
tion of Zc to eigenstates with with quantum numbers L, Lc and energy gE yields eigenstates 
with quantum numbers L -fl, + 1, gE. Similarly, the application of D^. to such eigenstates 
yields eigenstates with quantum numbers L — 1, L^ — l, gE. These properties are well known 
for the case of a two-body contact interaction |ll9| , |20| and generalize to the Hamiltonian (^. 
Note that the operators Ao and Ai are simply given by 

Ao = ^N{N-l), 

A^ = N{L-U). (8) 

Applying the operators Am onto the elementary symmetric polynomials yields 

Amei = Q for m > 1 (9) 

Thus, the one-dimensional space Vq = spanlef} is an invariant subspace of the Hamiltonian 
and therefore an eigenstate. Further, 

Ca = for m > 3. (10) 

This immediately shows that the (L — l)-dimensional vector space 

Vi =span{eL_Aet : A = 0, 1, 2, . . . , L - 2} (11) 

is annihilated by Am for m > 3. The space Wi = Vo U Vi is an invariant subspace of the 
Hamiltonian (^ since 

i2eA = 2AiVeA-2(iV-A + l)eieA-i, (12) 

and 

ii ca = ATA ca - (AT - A + 1) ei ca^i. 
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Eigenfunctions of the Hamiltonian (J^) in the subspace Wi are obtained by constructing 
states with good quantum numbers. To this purpose we project ei onto the space with 
= and obtain 



^ 1 



zn) 



m=0 



(13) 



The r.h.s. of the first hne is the Taylor expansion yielding a shift in the single-particle 
coordinates by —Zc- The second hne is obtained by combining the operators Dc and Zc to 
the operator for the angular momentum of the center of mass (El). This yields an operator 



that projects onto the subspace with Lc = 0. The expansion [117 



Zn - Zc 



L 

E 

m=0 



L - m , 



{-Z^r "'em{zu...,ZN). 



(14) 



shows explicitly that only states of Wi are involved in the construction of a state with 
Lr = 0. The states 



Zc ^L~\{Zi — Zc, ■ ■ ■ , Z]\i — Zc), 



X 



0, 1,2, ...,L-2,L 



(15) 



thus constitute an orthogonal basis of Wi. Since they have different quantum numbers 
Lc = X they are also eigenstates of the Hamiltonian (|^). Note that these wave functions 
do not depend on the coefficients Cm in the Hamiltonian. The corresponding energies are 
obtained by applying the Hamiltonian (|^) to the states (p!5|). Using eqs. (|),(|12]) and (p^lf ) 
yields 
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-N{N - l)co + N{L - A)(ci + 2c2 



A = 0,1,2,...,L-2,L. 



(16) 



The states (|T5|) thus form a ladder. Let us consider the important case of a two-body contact 
interaction. Analytical arguments showed that the ground state is obtained for X = L 
in the case of attractive interactions g < 0. For repulsive interactions, a combination of 
analytical |15|-0 and numerical calculations |]ll| showed that the ground state is obtained 
for A = 0. It is interesting to note that the corresponding fermionic problem, i.e. the 
quantum Hall effect with zero range interaction potentials, also has an exact ground state 



solution - the Laughlin wave function [21 



Note that the annihilation properties described in eq. (^) and eq. (^U[) may be generalized 
to further subspaces. Let 



span < CAi CAa 



ex, eT 



m 



0,1, 



L - 2k: 



k 



L 



m; 



A, >1 (17) 



denote the space of wave functions that contain products of k elementary polynomials dif- 
ferent from ei and cq. Using eq. @ and eq. ([T0|) one obtains 







for m > 2k + 1. 



(18) 
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However, these spaces are not invariant subspaces of the Hamiltonian, since e.g. 
62 (-21, ... , Zq) G Vel{zi, . . . ,zq). Let us consider the space that is generated when acting 
with the Hamiltonian onto Vk for k > 1. A direct calculation shows that 

{di - d^r i^m - Zr^nOm - dnY Vfe = for ^ > z/ + 2A; + 1. 

Therefore, 

AmVVk = ^ for m > 4fc + 1. 



Introducing the subspaces 



yields 



VWkdWok forA;>l. 



^ yv2k 

The Hamiltonian matrix thus can be cast into a block-tridiagonal form. The invariant 
subspaces Wq, and Wi yield a block- diagonal structure, and the remaining matrix is block- 
tridiagonal. 

The knowledge of L eigenstates (|15|) permits us to determine two additional eigenstates 
that are independent of the specific coefficients in the Hamiltonian (|^). For L = 4 and 
L = 5, there are two basis states with = 0. One of them has already been determined 
above, i.e. ei^Zi — Zc, . . . , zj^ — Zc). The other two wave functions in question are 

N / N 

N N N 

4>5 = NY^izj - Zcf + 2 Y.{zj - z,f Y.(zi - z^ 
j=l j=l 1=1 



i3 



These wave functions are polynomials of degree four and five, respectively. Therefore, one 
finds 04, 05 e VV2. 

The results derived so far were obtained for < L < A^. For L > N the Hilbert space de- 
pends on both L and N. The definition of the spaces Vk has to be modified correspondingly, 
e.g. Vi = span{eL_A : A = max {0, L — N), . . . , L — 2} to be compared with eq. 
However, the annihilation property described in eq. (|1^) remains valid and the Hamiltonian 
matrix continues to be of block-tridiagonal form. Note that the wave functions (|I5|) are 
restricted to parameter values A = max (0, L — A^), . . . , L — 3, L — 2, L for L > N. 

Quantum mechanical systems with invariant subspaces have received considerable atten- 
tion in recent years, since (quasi) exactly solvable systems fall into this class. For a general 
discussion see, e.g. refs. |p2|-p5|. In difference to the present work, the interest there is on 
m/inzte-dimensional Hilbert spaces with a nested sequence of finite dimensional subspaces, 
and most research has been limited to systems with a few degrees of freedom. Note finally 
that partially solvable many-body systems or systems with partial dynamical symmetry 



8l also permit the analytical computation of selected eigenstates. 
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In summary, we studied a large class of harmonically trapped, interacting boson systems 
at angular momentum L. The Hamiltonian leaves a L-dimensional subspacc invariant, and 
the quantization of the angular momentum for the center of mass yields several eigenstates 
that are independent of the details of the two-body interaction. Important examples of suit- 
able interactions are two-body potentials of zero range and interactions that are analytically 
in the squared inter-particle distance. 

This research was supported by the U.S. Department of Energy under Grant DE-FG03- 
00-ER41132. 
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